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ABSTRACT 

DR21(OH) is a pc-scale massive, -7000 M© clump hosting three massive dense cores (MDCs) at an early 
stage of their evolution. We present a high angular-resolution mosaic, covering ~70"xl00", with the IRAM 
PdBI at 3 mm to trace the dust continuum emission and the N2H + (J=l-0) and CH 3 CN (J=5-4) molecular 
emission. The cold, dense gas traced by the compact emission in N2H + is associated with the three MDCs 
and shows several velocity components towards each MDC. These velocity components reveal local shears in 
the velocity fields which are best interpreted as convergent flows. Moreover, we report the detection of weak 
extended emission from CH3CN at the position of the N2H + velocity shears. We propose that this extended 
CH3CN emission is tracing warm gas associated with the low-velocity shocks expected at the location of 
convergence of the flows where velocity shears are observed. This is the first detection of low-velocity shocks 
associated with small (sub-parsec) scale convergent flows which are proposed to be at the origin of the densest 
structures and of the formation of (high-mass) stars. In addition, we propose that MDCs may be active sites of 
star-formation for more than a crossing time as they continuously receive material from larger scale flows as 
suggested by the global picture of dynamical, gravity driven evolution of massive clumps which is favored by 
the present observations. 

Subject headings: ISM: kinematics and dynamics — stars: formation 



1. INTRODUCTION 

Two competing scenarios are challenged by observations 
to describe the formation of rich clusters and of high-mass 
stars: a quasi-static evolution (also known as core accretion 
model) versus a highly dynamical model. The first one is a 
turbulence regulated scenario, where a high level of micro- 
turbulence acts as an effective additio nal thermal pressure to 
balance gravity (e.g. McKee & Tan 2002). In the second 
one, massiv e cores form and evolve via highly dynamical pro 
cesses (e.g. |Bonnell & Bate |200 6 ; Vazquez-Semadeni et al. 
20071 IHeitsch etal. 2008 1 Hennebelle et al. 2008 ; Klessen & 



Hennebelle ||2010|), where convergent flows driven by large- 



scale turbulence, gravity and Galactic motions create dense 
structures down to sma ll scales by shock-dissipation at their 
stagnation points (e.g. Fiel d et aL||2008| Vazquez-Semadeni 
let al. |201 

In the Cygnus X region a population of massive dense cores 
(hereafter MDCs) was revealed by a systematic dust contin- 
uum survey ( [Motte et al. ||2007| ). Th e most massive of them 
are the birth place of high-mass stars (Bontem ps et al. |20T0] ). 
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Among the MDCs more massive than 40 M all show star for- 
mation activity, leading to the suggestion that the formation of 
high-mass stars within MDC s is a fast process ( [Motte et al 
2007 ). |Csengeri et al. ( 201 1) indeed recognized the important 



role of dynamical processes inside the young IR-quiet MDCs 
with crossing times only slightly larger than the local free-fall 
times, confirming a fast evolution. The precise origin of the 
MDCs and of their prime properties (mass, size, spatial dis- 
tribution) is still uncertain as well as one needs to understand 
how a rich cluster can be formed from relativel y small MDCs 
(see d iscussion in Bontemps et al. 2010). Csen geri et al. | 
( 201 1 ) also proposed that the massive protostars in the MDCs 
are found at the location of velocity shears by the convergence 
of small scale flows. The existence of these convergent flows 
needs however confirmation by, for instance, t he direct detec 
tionof the associated low-velocity shocks (e.g. Klessen et al. 
[2005] ). 

The clump associated with DR21(OH) is the most massive 
1 pc-scale clump of the whole Cygnus X region. It is embed- 
ded in the DR21 filament which contains as much as 30 % of 
the t otal mass in dense gas i maged by Motte et al. (2007) (see 
also Schneide r et al. |[20T0 ), and contains three 0.1 pc-scale 
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MDCs: CygX-N44, CygX-N48, and CygX-N38. CygX-N44 



hosts several spot s of OH, H 2 Q maser emissio n (Plambeck 
& Menten |[T990| |Liechti & Walmsley |fT997] ) and at high 



angular-resoluti on two strong peaks of continuu m emission 



were identified ( [Woody 
(MMl and MM2nT^ 
tains a "hot-core" (T>1 
sion (Mangum et al 



vo strong p i 
pdy et al.] ' 
Fig lib. The 



19891 |Mangum et al. |[1991) 



e brightest peak, MMl, con- 
K) and shows cm continuum emis- 
1991) suggesting the recent birth of 



high-mass star(s). CygX-N48 and CygX-N38 have no mid-IR 



emission and host cold gas (T<40 K) (Mangu m et al. |1 992). 

With such a high mass, the DR21(OH) clump qualifies as 
the best target in nearby regions (at less than 3 kpc J] to probe 
the initial conditions of a rich and massive protocluster. It is 
expected to form a cluster as rich or even richer than the Orion 
Nebula Clustei^] We mapped DR21(OH) at high-angular res- 
olution at 3 mm continuum, and in the N2H + (J=l-0), and 
CH 3 CN (J=5-4) molecular lines with the PdBI to trace the 
cold, dense gas and warm gas with complex chemistry, re- 
spectively. 

2. OBSERVATION AND DATA REDUCTION 

A mosaic with 3 pointings was obtained with the IRAJVj^] 
Plateau de Bure Interferometer (PdBI) with 6 antennas in Bq 
and Cq configurations between December 2008 and March 
2009. Baselines range from 24 to 452 meters and T sys was 
-100-150 and 60-80 Kelvins for the B and C tracks, respec- 
tively. The signal was correlated at 93.17613 and 91.98705 
GHz in order to obtain the N 2 H + (J=l-0) and CH 3 CN (J=5- 
4, K=0-4) lines respectively, with a velocity resolution of 
-0.25 km s" 1 . Line-free continuum emission was obtained 
with 6 broad band correlator units. We used as phase calibra- 
tor the bright nearby quasar 2013+370 and as flux calibrator 
the bright evolved star MWC 349. 

We used the GILDAS softwarf] for the data reduction 
and analysis and applied the same method for all datasets. 
Zero-spacing infor mation was available only for the N2H + 
(J=l-0) transition ( Schneide r et al. |20T0] ) and was combined 
with the interferometric data following standard procedures 
in GILDAS. ^ The image was restored with a natural beam 
weighting in order to favor sensitivity. Clean components 
were searched for within a polygon, which excluded the bor- 
ders of the mosaic with a higher noise level. The cleaning 
procedure was done with the Hogbom algorithm and compo- 
nents were searched down to 3<x per channel noise level. The 
resulting parameters are shown in Table [T] As a default in 
mosaic mode, the resulting clean maps are corrected for beam 
attenuation. 

3. N 2 H + AS TRACER OF THE COLD, YOUNG GAS 



N2H + is a tracer of cold, dense gas (Tafalla~et al. ||2002). 



Since it is quickly destroyed in the heated envelopes of proto- 
stars, it shall represent the starless/unprocessed gas, which has 
either not yet fragmented or originates from recently formed 
pre/protostellar cores. 



We adopt a distance of 1 .7 kpc (Schneider et al. 



2006). 
^Tthe 



2 Assuming a star formation efficiency of 30 %, the -7000 M Q of 
DR21(OH) would produce a stellar mass of 2100 M which is mor e than 
twice the stellar masses in the Orion Nebula Cluster ( Hillenbrand 1997). 

3 IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN 
(Spain). 

4 See http://www.iram.fr/IRAMFR/GILDAS/ 

5 We exploit here both the combined and the interferometric-only dataset 
of the N2H + data, because the latter one filters out extended emission from 
the clump and shows the compact emission, which is the main focus of this 
paper. 



In a quasi-static view of star formation, the pre-collapse 
core traced by N2H + is expected to evolve to a higher level 
of concentration by a progressive loss of support until it gets 
gravitationally unstable and collapses. In nearby low-mass 
star-forming regions such as p Ophiuchi, N2H + is indeed 
found to show single, narrow lines which are consistent with 
being close to or under gravitational collapse with only a 
small contribution from non-thermal motions ( [Andre et al. | 
[2007] ). 

3.1. Compact N^tt emission is concentrated in the MDCs 

In Fig. [I] a the spatial distribution of N 2 H + at low (IRAM 
30m) and nigh (PdBI) spatial resolution is compared with the 
dust continuum emission which is taken as a reference to trace 
the total column density. 

Seen with a single-dish telescope the distribution of N2H + 
differs from the dust continuum (Fig. [T^), as N2H + peaks in 
the south of the clump between CygX-N48 and CygX-N38. 
This prominent north- south inhomogeneity indicates that gen- 
erally the southern part of the clump is younger than the north- 
ern part. At high spatial resolution with the PdBI (Fig. [lj>), 
the compact emission of N2H + is clearly concentrated in trie 3 
known MDCs, CygX-N48, CygX-N38, and CygX-N44. This 
discrepancy of the distribution of N2H + seen with the single- 
dish telescope compared to the interferometer shows that a 
large fraction of the N2H + emission originates from extended 
structures that are filtered out by the interferometer. It also 
shows that the most compact components of N2H + , and there- 
fore the densest parts of this unprocessed, cold gas are asso- 
ciated with the MDCs in the clump. 

At high angular-resolution this spatial coincidence between 
the dust and N2H + and the fact that both tracers seem to form 
structures at similar size scales of ~0.1 pc further confirms 
the existence of a typical size-scale for MDCs and indicates 
that the MDCs are the main sites of present star-formation. 
In contrast, the hot-core region of DR21(OH)-MMl is de- 
void of N2H + suggesting that star formation there is in a later 
stage. A new, rather compact N2H + core is seen north of 
CygX-N38 (Fig.[lJ>). This new structure coincides with very 
weak (~1.3 mJy/beam peak intensity) continuum emission 
and could either be a fluctuation in the global diffuse N2H + 
gas, or could correspond to a new MDC in formation. 

3.2. Bulk motions in the N2H + gas associated with the MDCs 

At low angular-resolution with the IRAM 30m telescop e 
the spectra show a single component ( Schneid er et al. |2010] ), 
such as in p Ophiuchi, but with a large velocity dispersion of 
~1.3 km s" 1 . But while in p Ophiuchi th e line stays single at 
small scales ( |Di Francesco et al~||2004| ), the N2H + emission 
seen with the interferometer clearly splits here into individual 
velocity components (Fig. [2]). 

This is a striking feature and an important difference to 
low-mass star- forming regions. Since we could include the 
zero- spacing information from the IRAM 30m telescope, we 
are confident that there are no strong side-lobes and filtering 
which could modify the line profiles. Moreover, it seems that 
it is towards the MDCs that the individual velocity compo- 
nents can be recognized the best (see Fig. [2^). 

We show in Fig. [3] the maps of integrated intensity of the 
two velocity components which are associated with the two 
coldest MDCs, as well as position-velocity cuts through and 
across the intersection regions between the two velocity com- 
ponents. These cuts show velocity jumps indicated by arrows 
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Fig. 1. — a) Color scale shows the integrated intensity map of the N2H + (J=l-0) line obtained by Schneider et al. ( 2010) with the IRAM 30m telescope 
(integration range is between -10 to +0.5 km s _1 ). W hite contours indicate the 20% level of the peak increasing by 10%. Black contours show 1.2 mm emission 
obtained by the MAMBO-2 survey of Cygnus-X by IMotte et al. | f2007 1. Stars indicate sources with compact free-free emission (VLA archival data). Red 
crosses show DR 21(OH)-MMl, -MM2 and the peak 3mm PdBI continuum emission towards CygX-N38 and -N48, respectively. Labels mark the MDCs from 
Motte et al. (2007). b) Map of the same line obtained with the PdBI (without zero-spacings) integrated over the more isolated, F= 101-012 hyperfine component 
corresponding to -10 to +5 km s _1 , white contours indicate the 20% level of the peak increasing by 10%. Black contours show the 3mm dust continuum, c) 
Map of the 3mm dust continuum emission. Gray contours go from 5cr to 30<x by 5<x (<x = 0.3 mJy/beam), light grey contours go from 30cr to 150cr by 15cr. 
d) Map of integrated intensity of CH3CN emission over all K-components from -10 to +55 km s _1 . Gray contours go from 3<x to 20<x by 3<x, then by 40<x 
(<r = 35 mJy/beam). (For a color figure see the online edition.) 



which have relative velocity difference of 2 - 3 km s" 1 and 
are referred below as velocity shears. We note that the inter- 
section regions of the velocity shears are close to the strongest 
peaks of continuum emission, but do not coincide exactly with 
them. The velocity shears are found offseted from the strong 
continuum sources. These N2H + velocity patterns are very 
similar to what we have recentl y obtained using H 13 CQ + in 
isolated MDCs in Cygnus X ( [Csengeri et al. 2011| ). Simi- 
larly, we propose that the individual velocity components seen 
in N2H + are also best interpreted as flows converging to the 
gravitational well of MDCs. These flows may carry a sig- 
nificant amount of angular momentum which could lead to 
rotational motions. For a 200 M© core the Keplerian velocity 
at 0.1 pc is -2.9 km s" 1 which is at a similar order as the gra- 
dient observed in velocity field. On the other hand the pv-cuts 
in Fig. Pfl show that the maximum velocities are not found at 
large offsets from the velocity jumps similarly to what is ex- 
pected in a Keplerian rotation, but they show a more complex 
distribution. The observed velocity pattern of the gas here is 
thus not compatible with a pure rotation of the cores. Further- 
more an homogeneous sphere in rotation would also not show 
separated individual velocity components. 

It is surprising to see these flows in N2H + , since it normally 
traces only pre-stellar cores which are gravitationally bound 
structures. On the other hand, it may just indicate that the pre- 
stellar gas is already dense and cold in the convergent flows 
which then form the even higher density structures at the cen- 
ter of the MDCs. Furthermore we note that the observed range 
of velocity (between -6 .5 and +0.5 km s" 1 ov er the whole 
clump) is larger than in Csenge ri et al. ( 2011| ) most proba- 
bly because of the larger mass and gravitational well of the 
DR21(OH) clump. 



the local Jeans-mas^] which should concentrate in the cen- 
tral regions close to the gravitational center of the clump. In- 
stead, the distribution of dense gas traced by N2H + mostly 
reveals three centers of collapse, the MDCs. It is then not 
clear why the large scale, gravitationally driven flows split to 
predominately form a few MDCs dispersed over the clump. 
On the other hand, the observed hierarchical fragmentation 
is v ery similar to what is o btained in numerical simulations 
(e.g.fHennebel le et al. |201 1) . These models indicate that frag- 
mentation Ti^rovoCeaDya complex interplay between grav- 
ity, rotation, thermal/turbulent p ressure and magnetic forces 
( Vazquez-Semadeni et al. 2007 ) which may make it difficult 
to clearly pinpoint the physical origin of the observed frag- 
mentation properties. 

3.4. Are MDCs living for more than a crossing time? 

Since the gas dynamics dominate the evolution of the 
MDCs, the crossing time^Jof the MDCs, ranging from 5 to 
7 x 10 4 yr, should measure the life-time of the presently ob- 
served dens e gas. On the other han d the large scale flows 
observed by Schneider et al. (2010) are massive enough to 
continuously replenish the mass of the MDCs keeping them 
'active' for a longer period if the local gravitational wells are 
spatially stable. The presence of still active N2H + convergent 
flows at the location of each MDC favors this view of stable 
gravitational wells over time and thus suggests that the MDCs 
are re-filled (by large-scale flows) over longer times. 

It is also striking to see that all three MDCs are situated 
close to embedded IR-bright sources likely excited by em- 
bedded massive stars (see Fig. BJ. This reveals the presence 
of a population of young stars which could have been formed 
at earlier time inside the same three MDCs. 



3.3. On the origin of the clump fragmentation into MDCs 

If the -7000 M DR21(OH) clump would be governed 
only by thermal motions and gravity, it would fragment into 
a large population of ~0.4 M© fragments corresponding to 



6 The Jea ns-mass is calculated here for n ~ 1.3x1 5 cm 3 , and T = 
20 K {Motte et al. |2007) 

7 Since the crossing-times are almost equal to the local free-fall times 
which are ~ 5 x 10 4 yr, it re-inforces the view by which self-gravity is the 
main driver of the observed dynamics (the whole DR21 filament is in global 
collapse). 
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Fig. 2. — a) The map shows the 3mm continuum emission obtained with the PdBI (shown also in Fig.FTJ. The overlayed spectra correspond to the isolated, F = 
101-012 hyperfine component of the N2H + (J=l-0) emission (PdBI data combined with zero-spacings Trom the IRAM 30m telescope) resampled on a 4" grid. 
The range of the x velocity axis is from -8 to +0.5 km s _1 . The intensity range, y axis, scales from -0. 1 to 0.4 Jy/beam. Dots mark the position of a grid sampled 
on half of the beam. Green crosses point the positions where spectra are extracted and shown in panel b and c. b) and c) Extracted spectra towards the 3 positions 
in CygX-N48 and CygX-N38. Red and blue shading indicates the integration ranges shown in Fig. [3] (For a color figure see the online edition.) 



4. CH 3 CN REVEALS EXTENDED WARM GAS IN EACH MDCS 

4.1. Intriguing extended CH 3 CN emission in all MDCs 

A map of integrated CH 3 CN emission over all K=0-4 com- 
ponents is shown in Fig.[TJd. CH3CN has long been consid- 
ered to be produced only in hot cores an d hot corinos (e.g. 
Kurtz et al. 2000; Bottinell i et al. ||2004| ) as a second gener- 
ation molecule formed in the gas phase after ice sublimation 
resulting in very compact emission towards protostars. As 
expected, a bright, compact peak of CH 3 CN emission is cen- 
tered on the hot core source MM1 (Fig.[T]d). Among the other 
continuum sources, we do not detect any similar strong and 
compact CH3CN emission which would indicate hot cores. 

In addition to the compact emission of MM1, widespread 
extended emission of CH3CN is detected towards all MDCs of 
the region, including the colder ones, CygX-N38, CygX-N48. 
This extended emission is associated with the MDCs but is in 
detail almost entirely anti-correlated with the dust continuum 
(Fig.[3]b and e). In CygX-N48, it forms an elongated structure 
almost perpendicular to the dust emission. In CygX-N38, it 
is shifted by 5" from the dust continuum. This indicates that 
CH 3 CN has to be over-abundant in the gas (in comparison 
with the bulk of gas traced by the continuum) to explain this 
particular spatial distribution. 

4.2. Is CH3CN tracing shocks of the convergent flows? 

In Fig. |3]b and e, the two main N2H + velocity components 
are overlaid on the CH 3 CN emission. CH 3 CN seems to be 
bright at the location of the velocity shears (interfaces be- 
twee n th e red and blue velocity components) discussed in 
Sect. 3.2 and which we interpret as tracing N2H + convergent 
flows. It is striking that CH3CN coincides much better with 
these velocity shears than with the dust continuum. Therefore 
we suggest that CH 3 CN could trace warm gas associated with 
the shocks expected at the convergence of the flows. 



CH 3 CN has actually been detecte d towards a young pro - 
tostellar outflow shock in LI 157 ( [Codella et al. 1 12009). 
The abundance of CH3CN can have significa nt con tribution 
from grain surface chemistry (Garr od et al. ||2008] ) produc- 
ing CH 3 CN directly on ices. A higher CH 3 CN abundance 
could originate in warm gas from ice sublimation. The in- 
volved velocities are however here smaller than in outflows, 2 
to 3 km s" 1 in projection which converts to 3.7 to 5.5 km s" 1 
if corrected for an average p rojection anglaj From Fig. 3 in 



Kaufman & Neufeld (1996) we see that the resulting post- 



shock gas temperature could be close to 100 K for a magne- 
tized C-shock. At this gas temperature, the gra ins stay how- 
ever cold and ice sublimation is not expected (Dr aine et al. 
1983| ). Partial desorption has however been recently pro posed 
to explain gas phase abundance of complex molecules ( |Gar ' 



rod et al. ||2007| |Arce et ai~l|2008t [pberg et al. ||2009}. The 



higher gas temperature and sputtering of grain mantle in the 
shock may have led to an increased desorption of CH3CN 
which could explain the detectable abundance of CH3CN in 
the gas phase which is otherwise virtually not present in the 
cold, dense gas. 

5. CONCLUSIONS 

It is striking that N2H + , a classical tracer of cold, dense gas, 
shows flows associated with the highest density regions form- 
ing new high-mass stars. These flows are showing velocity 
shears which are found to spatially coincide with intriguing 
extended CH3CN emission. 

The evolution of the MDCs is driven by supersonic flows 
which are able to continuously supply a replenishment of ma- 
terial thanks to large scale flows (see Schnei der et al. | 2010) 
while the MDCs form stars. As a consequence, they mam- 

8 The true post-shock velocities could be higher if the flows accelerate in 
the gravitational potential of the MDCs. 
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Fig. 3. — a) Zoom into the CygX-N48 MDC. The map shows the 3mm PdBI continuum emission. Contours display the N2H + emission integrated between 
-1.1 to -3.9 km s _1 (red) and -3.9 to -6.9 km s _1 (blue) and go from 10cr increasing by 10<r. b) same as in (a) with the map showing the integrated emission 
over all K-components of CH3CN. c) Position- velocity cuts of the N2H + emission across (solid line) and along the shears (dashed line), d-f) Same as (a-c) for 
the CygX-N38 MDC with N2H + emission integrated between -3.9 to -0.9 km s _1 (red) and -6.7 to -3.9 km s _1 (blue) for the contours. 



tain active star-forming processes for longer than their free- 
fall times and crossing times if the convergent point of flows 
is stable in space over time. 

The spatial coincidence between the N2H + velocity shears 
and the intriguing extended emission in CH3CN indicates 
that we have probably detected for the first time the effect 
of low-velocity, but high density shocks expected in the case 



of widespread convergent flows building up new high density 
gas in the MDCs. 
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APPENDIX 



TABLE 1 

Parameters of the observations. For the molecular lines the rms noise is given for 0.25 km s _1 spectral resolution. 



PdBI 


synthesized beam 


P.A. 


rms 




["x"] 




[mJy/beam] 


3 mm continuum 


2.11 x 1.9 


84° 


0.3 


N 2 H + (J=l-0) 


2.09 x 1.89 


84° 


8.9 


CH 3 CN (J=5-4, K=0-4) 


2.13 x 1.92 


83° 


9.1 




Fig. 4. — Color composite image of Spitzer IRAC 3.6 jum, 4.5 jim and MIPS 24 /urn ( |Hora et al. |2009) . Contours show the 3mm continuum emission obtained 
with the PdBI. Note the bright embedded sources close to the MDCs. Note also that there is a faint 24 /urn emission at CygX-N38 which is an artifact of the MIPS 
PSF due to the nearby bright source. (For a color figure see the online edition.) 



